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TYPE II MIGRATION OF PLANETS ON ECCENTRIC ORBITS 
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ABSTRACT 

The observed extrasolar planets possess both large masses (with a median Msini of 1.65 Mj) and 
a wide range in orbital eccentricity (0 < e < 0.94). As planets are thought to form in circumstellar 
disks, one important question in planet formation is determining whether, and to what degree, a 
gaseous disk affects an eccentric planet's orbit. Recent studies have probed the interaction between a 
disk and a terrestrial planet on an eccentric orbit, and the interaction between a disk and a gas giant 
on a nearly circular orbit, but little is known about the interaction between a disk and an eccentric gas 
giant. Such a scenario could arise due to scattering while the disk is still present, or perhaps through 
planet formation via gravitational instability. We fill this gap with simulations of eccentric, massive 
(gap- forming) planets in disks using the hydrodynamical code FARGO. Although the long-term orbital 
evolution of the planet depends on disk properties, including the boundary conditions used, the disk 
always acts to damp eccentricity when the planet is released into the disk. This eccentricity damping 
takes place on a timescale of 40 years, 15 times faster than the migration timescale. 
Subject headings: planetary systems: formation 



1. INTRODUCTION 

The detection of over 200 extrasolar planets has over- 
turned our understanding of what constitutes a typi- 
cal planetary system. In contrast to the solar system, 
the extrasolar planetary systems generally contain giant 
planets with small semi-major axes (0.017 AU < a < 
6 AU) and sizabl e eccentricities (0 < e < 0.94) (e.g., 
IButler et al]|200l . 

The small semi-major axes of giant planets suggest 
that the many cxoplanets did not form in situ. Rather, 
planets are assumed to form exterior to the "snow-line," 
or ice condensation radius, and to subsequently migrate 
inwar d through inter actions with a circumstellar disk 
(e.g.. Ilda fe Linir2004h . This migration has two limiting 
cases in a laminar disk. In Type I migration, a planet 
lacks sufficient mass to clear a gap in the disk mate- 
rial and is driven inward, and i ts eccentricity damped, 
by a density wake in th e disk (j Papaloi zoii fc Larwoodl 
20001: IArtvmowic3 119931 : IPapaloizou fc Larwoodl 120001 : 
Cresswell et al.l 120071 ). In Type II migration, massive 



planets clear a gap in an annulus surrounding their or- 
bits. Analytic treatments of Type II migration gener- 
ally consider the planet to be locked t o the the viscous 
evolutio n of the disk ( e.g., lWardlll997[ ). although recent 
work by lEdgarl (|2008j ) shows that the planet is not as 
strongly coupled to the disk viscosity as predicted. He 
claims that the inverse relationship between migration 
rate and planet mass, combined with a weak dependence 
on disk viscosity, indicates that Type II migration should 
be thought of as an angular momentum exchange be- 
tween disk and planet, rather than the outer and inner 
disk using the planet as a simple relay station. 

The tendency of the planet to migrate inwards is well 
established, provided the plane t does not lie exterior to 
the majority of the disk (e.g. iLin fc Papaloizoul 119861 : 
IWardl [19971: lEdgarl [200l . The behavior of orbital ec- 
centricity under the influence of a disk, on the other 
hand, is not well understood and oft debated. While 



the aforementioned papers find that the disk damps 
eccentricity, some recent paper s argue that the disk 
can, in fact, excite eccentricity (IGoldreich fc Saril 120031 : 
IQgilvie fc Lubo^|[2003l : ID'Angelo et al.H2006D . However, 
previous studies started the planet on a circular or nearly 
circular orbit. Most known cxoplanets have substantial 
orbital eccentricity, and planets in multiple-planet sys- 
tems may attain significant eccentricities due to planet- 
planet interactions while Type II migration is still un- 
derway. Therefore, we turn our attention to Type II 
migration in the eccentric regime. 

Previously, iMoorhead fc Adamsl (|2008f ) performed cal- 
culations of Type II migration in the non-negligible ec- 
centricity regime and concluded that orbital eccentrici- 
ties up to values of 0.5 tend to be damped by the disk. 
In this paper we perform hydrodynamical simulations of 
a giant planet on an eccentric orbit. We describe our 
methods in iJ5J In tj3.ll we perform several tests to es- 
tablish the sensitivity of our results to simulation choices, 
such as resolution, and disk properties, such as viscosity. 
We present results showing that eccentricity is damped 
for all choices of disk properties in 33.21 and discuss the 
implications of these results in iJU 

2. METHODS 

All simulations presented in this paper show the orbital 
evolution of a Jupiter-mass planet in a minimum mass 
solar nebula. To perform o ur hydrodynam ic simulations, 
we use the FARGO code of lMassetl (|2000f ). which is a 2D 
polar grid code tailored to planet-disk systems. 

All simulations are performed with the standard pub- 
lic release of FARGO 1 , in which the disk is non-self- 
gravitating and isothermal. The hydrodynamics are 
largely c alculated in a "standard " manner (i.e., consis- 
tent with lStone fc Normari [19921 ). However, after calcu- 
lation of the radial transport and the residual azimuthal 
transport, FARGO separates the disk into rings and cal- 

1 http://fargo.in2p3.fr/ 
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culates the azimuthal transport in the comoving frame 
of each ring, before recombining the rings in prepara- 
tion for the next time step. This allows the code to take 
time steps larger than the Courant-Friedrichs-Lewy con- 
dition, which requires that disk material cannot move 
further than one grid cell per tim e step. A more detailed 
description of FARGO is given bvlMassetl d2000h . and an 
excellent synopsis is given bv lEdgarl (|2008| )~ 

Where sensibl e, we choose disk param eters equal 
to those in the Ide Val-Borro et all (|2006l ) comparison 
project. However, we must deviate from these choices 
in many cases where either: [1] the simulation of an ec- 
centric planet requires it, or [2] we wish to investigate 
the possible effects a disk parameter has on the long 
term evolution of the planet's orbit. As an example, we 
generally place our disk boundaries at 0.2 and 5.0 AU, 
rather than 0.4 and 2.5 AU, in order to accommodate the 
greater range in radial position of an eccentric planet. 

This work takes a n approach similar to that of 
iCresswell et al.l (|2007|) . in which the authors simulated 
the interaction between a disk and planets lOM^ in size 
with substantial initial eccentricity. As their work probes 
Type I migration, we focus on simulations with more 
massive planets (> Mj). Initially, we place the planet 
in a smooth disk (i.e., no gap), but integrate the system 
with the planet's orbital elements held fixed until the 
disk settles and the planet clears a gap. 

In all cases, our planet's mass is one Jupiter mass, and 
its initial semi-major axis is 1 AU. In each run the planet 
is held on a fixed orbit for a period of time (100 years, 
unless stated otherwise), and then released so that we 
may measure its response to the disk. This allows the 
disk to relax into a configuration that is consistent with 
the planet's eccentric orbit. This initial period with fixed 
orbital elements is reasonable when the semi-major axis 
damping induced by the disk takes place on a timescale 
longer than that on which eccent ricity is damped; eccen- 
tricity damping is calculated by iPapaloizou fc Larwoodl 
(2000J) to be a factor of 100 faster than semi-major axis 
damping. 

3. RESULTS 

We present two sets of results: First, we discuss the 
effects that varying disk properties has on planet-disk in- 
teractions and thus on the orbital evolution of the planet. 
Second, we identify common behavior over different disk 
types and discuss the ramifications for planet migration. 

3.1. Disk Properties 
3.1.1. Disk Viscosity 

First, we investigate the role of viscosity in planet-disk 
interactions using three choices: a constant viscosity of 
;/ = 10~ 5 (in units where a = G(M* + Alp) = 1, or 
v = 4.5 xlO 7 m 2 /s), and constant alpha- viscosities of a = 
10~ 3 and a = 10~ 2 . The torques resulting from these 
three simulations are shown in Fig. [TJ In each case, the 
total torque exerted by the disk on the planet fluctuates 
with diminishing amplitude over time. The timescale 
over which these fluctuations diminish depends on the 
disk viscosity: it is a few years for a = 10 -2 , a few 
hundred years for constant v = 10~ 5 , and greater than 
the simulation time of 500 years for a = 10~ 3 . 

For our subsequent simulations we choose a = 10 -2 . 
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Fig. 1. — Total torque exerted on the planet by the disk, for the 
cases of v = 10~ 5 (yellow), a = 10 -3 (blue), and a = 10 -2 (red) 
when a Jupiter-mass planet moves on a fixed, circular orbit within 
the disk. In future simulations we use a = 10 — 2 and let the disk 
relax for 100 years before releasing the planet. 

This choice is convenient, as the disk rapidly adjusts to 
the planet's presence in this case, but is also motivated 
by observations of T-Tauri disk s that are most accu- 
rately modeled with a = 10~ 2 (lHartmann et all 119981: 
Stepinskil 119981: iHersant et all 120011 : lLachaume et al l 
2003t I Andrews fc Williams! 120071 ) . This choice of a re- 
sults in fast migration; i.e., migration on a thousand-year 
timescale. 

3.1.2. Disk Boundary Conditions 

Next, we present the results of varying the boundary 
conditions between runs, while keeping the resolution 
constant. Extending boundary conditions from [0.4, 2.5] 
to [0.2, 5.0] to [0.04, 5.0] AU has a small effect on the 
values of the final semi-major axis and eccentricity of a 
planet starting in a near-circular orbit (Fig. [2]) but the 
overall orbital evolution does not qualitatively change. 
(Note that to maintain the same grid resolution, we in- 
crease the number of radial grid cells from 128 to 256 to 
384.) 

While a planet with low orbital eccentricity stays near 
a = 1 AU, we would like to study planets with apprecia- 
ble eccentricity; such a planet would pass through a range 
of radii of [0.7, 1.3] AU. Fig. [3] depicts simulations with 
the same set of boundary conditions as those shown in 
Fig. [21 the sole difference is that the planets of Fig. [3] be- 
gin with an eccentricity of 0.3. We find that when planets 
possess non-negligible eccentricity, the choice of bound- 
ary conditions significantly effects the long-term orbital 
evolution. Yet de/dt immediately following release of the 
planet is independent of grid size, and a(t) is similar for 
the two larger disk cases, suggesting convergence. 

The data presented in Figs. [T] through [3] are the re- 
sult of using FARGO 's "non-reflecting boundary condi- 
tions," which removes the reflected wake of the planet 
at the boundary, but does not allow the outflow of disk 
material. We perform analogous simulations using open 
boundary conditions where mass is allowed to flow off the 
edges of the grid. Again, we vary the position of these 
boundary conditions as depicted in Fig. [2l and again, 
we find that varying the boundary locations affects the 
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Fig. 2. — Semi-major axis and eccentricity evolution for a planet on an initially circular orbit at 1 AU with varying grid size, color coded 
to correspond with the diagram at bottom left. 



long-term orbital evolution, but not the initial de/dt, 
of the planet (Fig. 03) While these plots inform us as 
to the effects of varying the boundary conditions in hy- 
drodynamical simulations, we encountered a pronounced 
slowing-down of FARGO before these effects converged. 
With this limitation in mind, wc will continue to explore 
the effects of simulation choices using boundaries located 
at 0.2 and 5 AU and a radial grid resolution of 256. 

Animations of a subset of these simulations are avail- 
able online. 2 

3.1.3. Grid Resolution 

In this section wc demonstrate that our choice of grid 
resolution is adequate for calculating eccentricity evolu- 
tion; care must be taken to ensure that the system is 
not artificially circularized. In Fig. 0] we present a suite 
of six simulations of a Jupiter-mass planet orbiting for 
100 years on a fixed orbit with a semi-major axis of 1 
AU and orbital eccentricity of 0.3 which is subsequently 
released and allowed to migrate under the influence of 
the disk. We find that our choice of radial (N r = 256) 
and azimuthal (Ng = 384) grid resolution is adequate 
for our simulations; only when the number of azimuthal 
grid cells is cut to a third of this value do the results vary 

2 http:/ /www. astro. ufi.edu/~altheam/hydro. html 



significantly. All plots presented from this point forward 
will depict simulations with N r = 256 and Ng = 384. 

3.1.4. Surface Density 

We found that planet-disk interactions are directly pro- 
portional to the surface density, Co (where the surface 
density is given by a = ctq ( r / r o)~ p ), yet varying the 
slope of the surface density profile affects how long the 
torque oscillations take to settle down. The steeper the 
slope, the more persistent these oscillations are. This, in 
turn, alters the long-term orbital evolution of the planet's 
orbit, as can be seen in Fig. [5l For most disk profiles, 
the planet drifts inward upon release; however, in one 
case (p = —1/2) the planet initially experiences outward 
migration. In all other simulations we adopt p = for 
consistenc y with other hydrodynam ical studies, includ- 
ing that of Ide Val-Borro et all (|2006l ). 

3.1.5. Time of Release 

In previous sections we have reported results for sim- 
ulations in which we released the planet from a fixed, 
eccentric orbit at the time of apocenter. Here, we demon- 
strate that varying the mean anomaly at time of release 
does not affect the measured de/dt or da/dt. Fig. [S] shows 
the eccentricity as a function of time near the time of re- 
lease for each of these runs. The time of release in each 
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Fig. 3. — Semi-major axis and eccentricity evolution for a planet on an initially eccentric (eg = 0.3) orbit at 1 AU with varying grid size, 
color coded to correspond with the diagram in Fig. [2] Simulations shown in the top two plots have non-rcflccting boundary conditions, 
while simulations in the bottom two plots have open boundary conditions. 
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Fig. 4. — Orbital evolution for varying azimuthal and radial grid resolution. The disk boundaries are located at 0.2 and 5.0 AU. The 
planet is kept on a fixed orbit for 100 years, at which point the planet begins to migrate under the influence of the disk. The simulations 
are coded as follows: Red curves have a radial grid resolution of 256, blue curves have a radial grid resolution of 128. Simulations with an 
azimuthal grid resolution of 384 are solid, 256 are dashed, and 128 are dotted. 
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of these runs differs by one quarter of an orbital period 
and this difference produces small differences in e(t) and 
a(t). However, these differences arc indeed small; the 
right panels of Fig. [6] show that over tens or hundreds of 
orbits, the differences in orbital evolution resulting from 
varying the release time arc indistinguishable. 

Additionally, we can see from the left panels of Fig. [5] 
that while the semi-major axis and eccentricity appear 
to be smoothly damped in the lower two panels of Fig. [6l 
this is in fact not the case on sub-orbital timescales. 
When we calculate eccentricity and semi-major axis 
damping (or excitation) timescales for our simulations, 
we refer to the orbital evolution averaged over several 
orbits. 

3.2. Orbital evolution of an eccentric planet 

Our primary goal in this study is to determine the or- 
bital evolution of an eccentric planet in a thin, flat disk. 
In the previous section, we presented results of simula- 
tions with a variety of numerical parameters and disk 
properties, focusing on the semi-major axis and eccen- 
tricity evolution. We found that eccentricity damping 
occurred upon release of the planet regardless of resolu- 
tion, disk size, time of release, and surface density profile. 
We also note that the initial value of de/dt appears to be 
relatively independent of disk properties and numerical 
parameters, and thus we may draw conclusions that are 
insensitive to these quantities. However, changing the 
disk's properties - in particular, the location and type of 
boundaries - can substantially alter the long term orbital 
evolution of an eccentric planet. 

In this section we calculate the initial damping rate 
(de/dt) as a function of initial eccentricity for an eccen- 
tric planet-disk system. Previous plots show one of two 
cases, e = or e = 0.3. In Fig. [3 we present de/dt as a 
function of eo- The simulations depicted have the follow- 
ing properties: a constant alpha- viscosity of a — 10~ 2 , a 
flat surface density profile (p = 0), non- reflecting bound- 
aries at 0.2 and 5.0 AU, and 256 radial and 384 azimuthal 
grid divisions. Wc find that eccentricity is damped on a 
fairly uniform timescale, r e « 40 years. This is longer 
than an orbital timescale, but shorter than the migra- 
tion timescale r a , justifying our approach of holding the 
planet on a fixed orbit for 100 years while the disk re- 
laxes. 

The eccentricity damping timescale (r e 40 years) is 
independent of initial eccentricity and roughly 15 times 
faster than the migration timescale of a planet in a cir- 
cular orbit; the instantaneous value of (de / dt) / (da / dt) 
upon release varies widely due to varying initial values of 
da/ dt] some simulations produce an initial outward drift 
before inward migration takes place (see, for instance, 
the p = —1/2 case in Fig. [5]) 

4. CONCLUSIONS 

This study considered hydrodynamic simulations of 
Jupiter-mass planets in a disk with varying initial orbital 
eccentricity. In our simulations, the planet was initially 
held on a fixed orbit for a period of 100 years, allowing 
a gap to form around the planet's orbit before measur- 
ing da/dt and de/dt. Such a scenario could arise for a 
single planet in a disk due to the following processes: 
[1] The planet my form via gravitational instability in 



an eccentric orbit (Boss 2000). [2] The planet may be 
scattered by another body before the disk dissipates and 
forming a gap in the disk before its eccentricity is com- 
plete l y damped dChatteriee et al.l 120081: iMoeckel et al.1 
l200l iThommes et al.ll2008HMatsumura et al.ll2009l ). [3] 
The planet may enter a mean motion resonance with 
another planet; pairs of planets migrating inward while 
locked in a mean motion resonance experience an in- 
crease in orbital eccentricity until an equilibrium is 
reached between the eccentricity pumping of the migrat- 
ing resonance and the eccentricity damping of the disk 
(|Lee fc Pealdl20f)l IKlev et aLll2005h . 

Wc find that the choice of disk parameters can strongly 
affect the long-term (i.e., > 10 yrs) orbital evolution of 
the planet. For instance, low disk viscosity produces long 
lived oscillations in the disk torque. The placement of the 
boundaries affects whether and when a resurgence in or- 
bital eccentricity occurs, although this resurgence is less 
pronounced for open boundary conditions than it is for 
non-reflecting boundary conditions. If a very low grid 
resolution is used, the planet's orbit is artificially circu- 
larized, but this is not a concern for the grid resolution of 
N r = 256, N g = 384 used for this study. Finally, we find 
that the shape of the surface density profile affects the 
amplitude and period of oscillations in eccentricity af- 
ter release, as well as the extent and timescale of orbital 
migration. 

Despite the variety in the orbital evolution of the 
planet, we find that the initial drop in orbital eccen- 
tricity (on timescales of a few to tens of orbits) is con- 
sistent across the tested disk parameters. Furthermore, 
the eccentricity decay timescale r e = — e/ '(de/dt) at the 
time of the planet's release is equal to 40 years and is 
independent of starting eccentricity for eccentricities up 
to 0.5 (Fig. [7J. Simulating planets with eccentricities 
greater than 0.5 was not possible with our chosen disk 
size. This value of 40 years is similar to the timescale 
for eccentricity damping in Type I migration obtained 
by iCresswell et all (|2007l ). We measured r a /r e (where 
r a = —a/ (da/dt)) to be < 15. While eccentricity damp- 
ing took place on a shorter timescale than orbital migra- 
tion, this falls short of the value of 40 requir e d to repro- 
duce the GJ 876 system (|Lee fc Peald I2002t IKlev et al.l 
2005). This scenario of rapid eccentricity damping and 
less rapid inward migration is not dissimilar to Type I mi- 
gration, and thus suggests no reason for a difference in 
orbital element distribution between gas giant and rocky 
planets. 

Most hydrodynamical treatments of planet-disk in- 
teractions during Type I and Type II migration result 
in a net decrease in eccentricity, and no hydrodynam- 
ical treatment increases the plan et's orbital eccentric- 
ity beyond 0.15 (see, for example, iD'Angelo et al.l [20061 : 
ICresswell et aD [2003). Thus, we conclude that disk- 
planet interactions alone arc unlikely to excite the ec- 
centricities of exoplanets and to produce the observed 
range in eccentricity for giant planets. 

Planet disk interactions may be important for damp- 
ing eccentricities during the formation of multiple- 
planet systems; eccentricity damping is often re- 
quired to reproduce the observed range in orbital ele- 
ments dMoorhead fc Adamsll2005tlChatteriee et~aT1 l2008l: 
IThommes et al.l 120081: iMatsumura et al.l I2009il Fur- 
thermore, planet-disk interactions may be important 
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Fig. 5. — Orbital evolution for varying p. where <r = ao(r /ro)~ p ■ The disk boundaries are located at 0.2 and 5.0 AU. The planet is kept 
on a fixed orbit for 100 years, at which point the planet begins to migrate under the influence of the disk. The curves are coded as follows: 
red: p = —1/2, orange: p = 0, yellow: p = 1/2, green: p = 1, blue: p = 3/2. 
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Fig. 6. — Semi-major axis (left) and eccentricity (right) evolution of a Jupiter mass planet which is released at different phases in its 
orbit. Plots on the top show orbital evolution near the time of release, while plots on the bottom show orbital evolution for the full span 
of the simulations. Until 99.3 years have elapsed, the four runs are identical; the runs feature planets released from a fixed orbit at times 
separated by one quarter of an orbital period. 
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if these processes can reproduce the observed distribu- 
tion in cxoplanct orbital elements. Additionally, we 
will investigate the degree to which self-gravity affects 
planct-disk interactions in the eccentric regime, using a 
recently ava ilable version of FARGO that includes disk 
self-gravity (jBaruteau fc MasseHl2008D . 
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Fig. 7. — Eccentricity damping timcscale as a function of starting 
eccentricity for a set of hydrodynamical simulations (dots). The 
average value of this timescale for all shown simulations is underlaid 
as a dashed line. 

for setting the fin al orbital configuration of resonant 
planetary systems dLee fc Peald 12001 Kiev et all 120051: 
ISandor fc Klevl 120061 : ISandor et all 120071) . Future re- 
search should consider a combination of planct-planct 
interactions and planet-disk interactions, using our pre- 
diction damping timescales of r a /r e = 15, to determine 
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